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1 Introduction 



This article completes our investigations of the finite-size scaling properties of reaction- 
diffusion systems. Before getting into details we would like to outline a couple of questions 
which we will discuss in the present work: 

• Universality of the scaling properties 

To what extend are the scaling exponents and functions and the corrections influenced 
by special values of the reaction rates? 

• Influence of the initial conditions 

Does the choice of the initial condition change the scaling exponents or functions or 
does it enter only in the corrections? 

• Numerical methods 

Is it practically possible to compute the scaling properties by Monte Carlo simulations? 
Can extrapolations of exact finite lattice data provide information about the scaling ex- 
ponents and functions or even about the corrections? Exact numerical diagonalisations 
of the quantum chain Hamiltonians corresponding to a certain model are only possible 
for very small lattices. Typically one can handle about 20 sites for a two state and 10 
sites for a three state model. Therefore a point of major interest is to figure out what 
sizes of lattices are necessary to obtain reliable results. Furthermore we discuss which 
parameters such as boundary conditions, initial conditions, reaction rates or paricular 
values of the scaling variable z are most suitable for numerical investigations. 

We choose the single species coagulation model as an example to study these questions. 
In the previous two papers (hereafter referred to as I and II) we used analytical methods 
to investigate this model with periodic [||] and with open boundary conditions and for a 
very special value of the coagulation rate. We have shown that in the finite-size scaling limit 
L — > oo, t — > oo with the scaling variable z = kept constant the concentration of particles 
takes the form: 

c(z, L) = L x [F (z) + L-yF(z) + ••■]. (1.1) 

The scaling exponent x was found to be equal to —1 for both kinds of boundary condi- 
tions. The scaling functions F (z) however are found to be different for periodic and for open 
boundary conditions. 

We briefly review the definition of the coagulation model (for details we refer to paper I). 
The model is defined on an one- dimensional lattice of length L. The following reactions are 
allowed: 

• diffusion with rate D 

A + <&< — >$ + A, (1.2) 

• coagulation with rate c 

A + A — ► A . (1.3) 

We always choose units of time so that D — 1. In |§ it was shown that the master equation 
describing a recation-diffusion system can be mapped onto an euclidean Schrodinger equation. 
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Then the time evolution of the system is determined by the Hamiltonian of an one-dimensional 
quantum chain: 

||P) = -H\P) . (1.4) 

In paper I we have shown that the coagulation model is described by the Hamiltonian H = 
H + Hi where: 

Ho = ~ y t[ofof ¥1 + <rM +1 + A'afat +1 + (1 - A')K + af +1 ) + A' - 2] (1.5) 

and 

Hi = ~(l- A') X>+(1 - + (1 - ^Kil • (1-6) 
z i=i 

By: 

A' = l-c (1.7) 

we denote the asymmetry between the diffusion and the coagulation rate. This notation has 
already been used in paper I and ||. The spectrum is determined by H alone, which is the 
Hamiltonian for the XXZ-chain in an external magnetic field. As long as no decoagulation 
reaction (A — > A + A) occurs the system is on the Pokrovsky-Talapov line independently 
of A'. 



The first aim of the present paper is to investigate the scaling behaviour (O) of the 
concentration for different choices of the coagulation rate c, i.e. different values of A'. It 
can be shown 0] that in the continuum, along the Pokrovsky-Talapov line, the system is 
always massless with the same quadratic dispersion relation. We expect that this fact is 
reflected in the scaling properties. For A' = the spectrum can be obtained in terms of free 
fermions. This fact allows the analytical solutions presented in paper I and II. For A' 7^ 
these calculation do not apply. Therefore we use Monte Carlo simulations to determine the 
scaling exponent and function as well as the corrections. 

Our second point of interest is to figure out in which way the choice of initial conditions 
influences the scaling properties. For A' = the case of uncorrelated initial conditions with 
arbitrary initial occupation probability p was treated analytically in paper I and II. Now we 
use Monte Carlo simulations and extrapolations from finite lattices to study initial states with 
small clusters. These studies are carried out for A' = only. 

We test the accuracy of Monte Carlo simulations in the case of A' = where analytical 
results are available. The comparison of Monte Carlo data and analytic results demonstrates 
that Monte Carlo simulations provide reliable information about the scaling exponents and 
functions as well as about the corrections. 

It has been pointed out in paper I that scaling relations allow to compute estimates 
for critical exponents of an infinite system from the scaling behaviour of a finite one. An 
important purpose of the present paper is therefore to figure out what sizes of lattices are 
necessary to compute the scaling properties numerically. The case of A' = is used to test 
the accuracy of finite lattice extrapolations. It is shown that finite-size analysis allows to 
determine the scaling properties already from numerical data for lattices of length smaller 
than ten sites. 
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We decided to organize the paper according to the different numerical methods we use. In 
Section || we discuss Monte Carlo simulations of reaction-diffusion systems and present the 
results obtained by this method. Section ^ is devoted to extrapolations from finite lattices. 
The accuracy of the extrapolations is studied and the investigation of different initial condi- 
tions in the case of open boundary conditions is presented. We close with the discussion of 
our results. 

2 Monte Carlo Simulations 

Before we present the results of the Monte Carlo simulations we briefly explain how the 
simulations are performed. The dynamics of the system are determined by the rates r Qj/ g(7, 5) 
for the nearest neighbour interactions: 

(a,P)— ► ( 7 ,5). (2.1) 

Here a value of a, /3, 7, 5 = 1 corresponds to a particle and a value of to a vacancy. (This 
notation is connected to the one used in paper I by w a -.- /t p-s( i y, 5) = r at p(j, 5), the change is 
only made for convenience.) 

We consider two types of initial conditions: in an uncorrelated initial state all sites are oc- 
cupied with equal probability p at t = 0; a weakly correlated state consists of an interpolating 
sequence of two particles and two vacancies (AA$$AA$$ ■ ■ ■). 

We only study vacuum driven processes, i.e. no reactions occur on a pair of empty sites. 
Therefore we use the "direct method" ||, 1 for our simulations. By At we denote a discrete 
time step. Each Monte Carlo step consists of the following operations: 

1. ) A particle is chosen at random and afterwards, with equal probability, one of its neigh- 

bours. The resulting pair (a, 0) cannot consist of two vacancies. 

2. ) A new configuration (7,5) 7^ (a, (3) is chosen with the probability: 

n{a,P)T a ^{j,8)At (2.2) 
or the pair remains unchanged with probability: 

l-n(a,(3)At £ T a ^,5) . (2.3) 

Here n(a,(3) is defined by: 

f 1 if (a, (3) = (1,1) 
n{a,/3)=\ 2 if (a, (3) = (1,0) . (2.4) 
( 2 if (a, (3) = (0,1) 

This factor |7| simply accounts for the fact that a pair containing one vacancy is chosen 
with half the probability of choosing one that consists of two particles. 

3. ) Finally the time is increased: 

t — ► t + At/N (2.5) 
where N is the total number of particles on the lattice. 
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If the coagulation process has been chosen in a Monte Carlo step, the average concentration 
decreases by an amount of: 

* = zk <2S) 

where N s denotes the number of samples (i.e. program runs) used in the simulation. Therefore 
the simulated time t should not reach values where the concentration c(t) becomes of the order 
of q c . Because we are interested in the scaling behaviour and the corrections, it is necessary 
to perform the simulations on small lattices (L < 100). Therefore the number of samples N s 
has to be large. A large number of samples is advantegous because the standard deviation of 
the mean decreases like (iV s )~ ' 5 . The limit is set by the growth of the necessary CPU time. 
Averaging over N s = 20, 000 samples provides very accurate results as the test described in 
the sequel will show. 

The time discretisation At clearly influences the quality of the simulations. For smaller 
lattices this influence grows, because the number of particles is small (cf. Eq. ( |2.5| )). As a 
test, we perform simulations for A' = and periodic boundary conditions and compare the 
results with the exact expressions from paper I. The difference of the concentration averaged 
over 20,000 program runs and the exact value is used as a measure for the accuracy of the 
simulation. Tab. [I] shows the maximum value of the relative deviation which is observed 
during simulations up to a value of z max = 0.1. Then the maximum simulated time is 
tmax = z maxL 2 /4:. It can be seen that At has to be decreased to obtain results with equal 
accuracy for smaller lattices. The table shows typical values for At which are used for different 
lattice lengths. 

However, for very small values of At the probability that "something" will happen in 
a Monte Carlo step becomes very small in comparison with the probability that "nothing" 
will happen. Therefore the fact that a random generator never provides perfectly uniformly 
distributed numbers limits the possibility to improve the simulations by decreasing At. 

Another serious limitation is the growth of the necessary CPU time which is proportional 
to (At) -1 . As an example take the CPU time of 80 seconds for L = 3, N s = 20,000, 
At = 5 • 10 -3 and z max = 0.1. 

We are interested not only in the scaling function but in the corrections as well. The latter 
can only be determined for values of the scaling variable z for which the corrections are much 
larger than the numerical errors. Values of z — 0.01, . . . , 0.1 turn out to be suitable. 

To demonstrate the accuracy we make simulations for L — 9, 10, . . . , 16 with At = 5 • 10~ 3 
and averaged over 20,000 samples. The difference between the simulated concentration and 
the exact expression is always less than 0.4% and less than the standard deviation of the 
mean. With the help of a x 2_ fit we compute approximations for the scaling function and the 
corrections in for fixed values of z. The results are given in Tab. [TJ along with the exact 
expressions from paper I. It can clearly be seen that Monte Carlo simulations provide reliable 
information about the scaling properties. 

The programs were written in Fortran using the RAN random number generator and 
executed on DEC work-stations. 
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2.1 The Influence of the Coagulation Rate 

To solve the coagulation model exactly in paper I and II it was necessary to choose the 
coagulation rate equal to the diffusion rate, i.e. A' = 0. Now we use Monte Carlo simulations 
to study the way the scaling properties change for different choices of the coagulation rate, 
i.e. A' ^ 0. 

We use a full lattice as initial configuration corresponding to an occupation probability 
of p — 1. Both, periodic and open boundary conditions are considered. Averages are always 
taken over N s = 20, 000 samples, except where other numbers are stated. 

The first step is to determine the scaling exponent x in Eq. ( |TTT|) . For this purpose, for 
each value of A', we perform simulations for L = 500, 550 . . . 950 taking At = 0.5. Averaging 
over 400 samples is sufficient in this case, whereas in all other simulations we use 20, 000 
samples. For each value of z, x is determined with the help of the \ 2 test, according to: 

log (c(z)) = const. + x log(L) . (2.7) 

With an accuracy of 5% or better we obtain x — — 1, independently of A'. 

The second step is to find the leading term of the corrections. In order to do so, we need 
an approximation for the scaling function Fq. For this purpose we use the value of Lc(z, A') 
obtained from simulations on a lattice of length L = 1000. Then we perform another set of 
simulations on lattices of length L = 30, 35 ... 85 with At = 5 ■ 10~ 3 . Using this second set of 
data, we make a logarithmic fit of the type 

log (L c{z) - F Q ) = const. - y log(L) . (2.8) 

In all cases the correction exponent y is found to be equal to +1 with an error of 15% or less. 
Observe the contrast to the case of A' = where the leading correction exponent is y — 2. 
So the dominant correction for A' ^ is of the order ofl/L. 
Therefore we finally approximate the second set of data by: 

lcM1 . FlMH «Ml + *^. (2 . 9) 

The scaling function Fq is found to be independent of A' and has the values given in paper 



I, respectively II, i.e. F (z, A') = F (z } 0). This can be seen in the Tabs, [ill] and |VI|, where 
the values obtained for F are given for periodic and for open boundary conditions. 
The first order correction function Fi behaves according to the following scheme: 

Fi(z, A') > for A' > 
Fi(z, A') = for A' = 

F 1 (z,A') < for A' < . (2.10) 

The values for the first and second correction functions are given in Tabs. [TV] and [V] 



for periodic boundaries and in Tabs. [VII] and [VII]] for open boundaries, together with the 
analytic results for A' = 0. As the second order correction term, i.e. F 2 {z 1 A')/L 2 1 becomes 
already of the order of the statistical errors of the simulations the correction function F 2 can 
only be determined with large numerical errors. The correction functions are presented in 
Figures 1-4. 
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As has been explained in paper I the scaling limit for small z and the long-time behaviour 
in the thermodynamic limit (L — > oo) are related. For periodic boundaries and A' = we 
found the following expansion for the concentration: 



c(t) 




1 




+ ■ 



(2.11) 



The leading term is obtained from the scaling function and the second term from the expansion 
of the corrections in 1/L 2 . An analysis of the leading corrections for A' ^ which are of the 
order of 1/L and behave like 1/z leads to an expansion: 



c(t) 




a(A'p) 

1 + 1 ' + 



(2.12) 



Thus the change of the coagulation rate does not affect the leading term of the large time 
expansion of the concentration in the thermodynamic limit. However the next-to- leading 
term is no longer of the order of t~^ as for A' = but of the order of t~ l . The coefficient 
a depends on the value of A' and on the initial condition. For p = 1 and periodic boundary 
conditions we present the values of a for different A' in Tab. [IX|. These data are computed 



using a x'-fit of the values for the first correction function F\ (cf. Tab. IV). 



If one is interested only in the time dependence for systems in the thermodynamical limit 
(L — > oo) there is of course no need to to calculate the scaling expansion first. Instead one can 
directly perform simulations on large lattices. For A' = 0.5 we check the result given in Tab. 
|PX| by a simulation on a lattice of 2000 sites. A x 2_ fit of the data averaged over N s = 20.000 
samples gives a(0.5, 1) = 0.818 ± 0.006. This value is compatible with the one obtained from 
the scaling data. However the errors for the direct simulations are smaller by one order of 
magnitude. 



2.2 The Influence of the Initial Conditions 

For periodic boundaries and A' = we use Monte Carlo simulations to study the influence 
of a weakly correlated initial state. The configuration at t = is chosen to be a interpolating 
sequence of two particles and two vacancies (AAMAAM . . .). 

Lattice lengths of L = 36, 40 ... 80 (At = 5 ■ 10~ 3 ) are used. We analyse the data using the 
procedure described in the previous section. The results, shown in Tab. [X|, confirm that the 
scaling function is not influenced by an initial state with small clusters. The corrections are 
still of the order , but they are different from those given in paper I for p = 0.5, although 
the initial concentration is the same, c(0) = 0.5. 



3 Extrapolation from Finite Lattices 

The coagulation-decoagulation model, being exactly solvable, provides a good experimen- 
tal field to test the accuracy of numerical extrapolation from finite lattices to determine scaling 
functions and exponents. Our special interest is to find out from what size of lattices it is 
possible to obtain reliable information about the scaling exponent and the scaling function of 
the particle concentration as well as about the corrections. 
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A common method used to perform a finite-size extrapolation was proposed by Bulirsch 
and Stoer || (in the sequel being referenced as BST). 

We first shortly review this algorithm, following ||. Let T(r) be a function with an 
expansion 

Til) = T+ (IY T' + (lY T" + ... (3.1) 

oo. A table of extrap- 



.LJ \LJ \L, 

where we are interested in the asymptotic behaviour in the limit L 
olants is formed 
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(3.2) 



by application of the following rules 
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-1 



. (3.3) 



Here u; is a free parameter. The algorithm arises by approximating the function T (j^j by a 



sequence of rational functions in the variable 

<L) 



1 \u 



t[ L ^ is expected to converge faster than 



T to the limit T(0) when L — > oo. The parameter is determined by the condition that 
the errors of the approximation, i.e. the difference of two successive extrapolants T^ +1 ^ — Tffl 
are minimized for each m. This method is very convenient for our purpose, because a high 
precision for the extrapolated functions is obtained using only a few input values (which, 
however, must be known very precisely). 



3.1 Extrapolations for Periodic Boundary Conditions 

As we are interested to know how good the extrapolation from small lattices are, we 
make calculations for the scaling function, the scaling exponent, the correction exponent 
and the correction function for the coagulation model with periodic boundary conditions 
and for A' = 0. For the extrapolations, we take lattice lengths between 2 and 9 sites. We 
tried to work with seven input values (L = 2,3,4,5,6,7,8) and with eight input values 
(L = 2, 3, 4, 5, 6, 7, 8, 9). For small values of z the algorithm is convergent only in the latter 
case. The input values are calculated using the exact expression for the concentration in the 
coagulation model derived in paper I (cf. eq. (7.2)). 

For z the values z = 0.1, 0.2, . . . , 1.0 were chosen, covering a sufficiently large range in time 
between t = 0.1 and t = 25 where the behaviour of the scaling function and its corrections 
are investigated. 

Let us now describe the extrapolations in detail: 
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a) Scaling exponent 

In order to get a higher precision in our estimates, we write the scaling exponent as 
x = —1 + x and try to determine the difference x to the theoretical value for the scaling 
exponent x = 1, using the following expression to form the extrapolants: 



ln( (L+1) < L+1 ^ ' 
lnO L 



4 L) = ■ (3-4) 



The only inconvenience of this method is the fact that one obtains one input value less 
than for the extrapolation of the scaling function. The results are shown in Tab. fKl] for 
each z. We find x = with a numerical error less than 3 • 10~ 3 . 

b) Scaling function 

The scaling function F — 1 is extrapolated from 

T (L) = L c(z } L) — 1 . (3.5) 

The results are compared to the exact values of the scaling function F (z) — 1 calculated 
from paper I in Tab. |X1]. The bigger the values of z are, the higher the precision of the 
extrapolants becomes: for small z, the difference to the theoretical value is of the order 
of 2 • 10~ 3 , for z = 1.0 we even get a precision of 10 digits in comparison to the exact 
scaling function. This can easily be understood looking at the definition of z = ^jjr- A 
higher value of z means that the time t is already much later and that we are closer 
to the finite-size scaling regime where t — > oo. So the values of L c(L, z) for bigger z 
should better approximate the theoretical scaling function. 

c) Correction exponent 

Similar to a), we extrapolate the difference y = 2 + y to the theoretical correction 
exponent y = —2. The input values are calculated according to 



ln( (L+1) J (z,+ 1 "' / - r l - :| - r " ! - 1 



~(L) _ X11 V L*[L c{L,z)-F (z)\ ,„ R , 



For the value of the scaling function Fq(z) we use the result obtained in b). So there 
is already a numerical error in the input values leading to a lower precision for the 
correction exponent in comparison to the scaling exponent. We find y = with an error 
of 10~ 3 . The extrapolated value at z — 0.6 is worse because the correction function 
F 2 (z) is close to its zero. 

d) Correction function 

We are interested in the corrections F 2 (z) to the scalingfunction, i.e. we extrapolate 

T (i) = L 2 (Lc(L,z)-F (z)) . (3.7) 

For F (z) we again use the result obtained in b). Arguing similarly to the approximation 
of the scaling function, we can understand the high precision of the extrapolated correc- 
tion function for big values of z. But even for small values of z (except for z — 0.1) the 
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precision is in the range of 6 • 10~~ 4 . The bad value for the extrapolated corrections for 
z — 0.1 can be understood looking at the corresponding extrapolated scaling function 
which is known only with a precision of 2 ■ 1CT 3 . As this value enters the input values for 
the extrapolation of the corrections, the algorithm cannot lead to a reasonable result. 

In conclusion, we can say that the BST algorithm provides a powerful tool for extrapo- 
lations from finite lattices. Already from ten sites, for which the Hamiltonian can easily be 
diagonalized explicitly in order to get a value for the concentration, we obtain extrapolations 
with a very high precision. This knowledge can be used especially to examine non-integrable 
chemical models and to determine their critical exponents and finite-size scaling properties. 

3.2 Extrapolations for Open Boundary Conditions 

Once having tested the accuracy of the BST-algorithm we now use it to examine the 
influence of a weakly correlated initial state on the scaling and the correction function for 
open boundary conditions. For this purpose, we numerically calculate the concentration for 
finite lattices using the expressions found in paper II for the eigenfunctions of the system of 
differential equations describing the time evolution of the coagulation model (cf. eqs. (2.13)- 
(2.16)). Analytically, these calculations become too cumbersome. Afterwards, we extrapolate 
the finite lattice data with BST. This time we use for the input values lattice lengths between 
20 and 40 sites in order to get a higher precision for the extrapolated values. 

Using the same procedure as described in the previous section, we determine the scaling 
function and the corrections for the weakly correlated initial condition A400A400 . . .. The 
results are shown in Tab. [X1J| in comparison to the analytical values for uncorrelated initial 
conditions (initial occupation probability p = 0.5). 

The scaling function is again found to be independent of the initial conditions with a 
precision of 6 digits for values of z larger than z = 0.1. 

On the other hand, the corrections are influenced by the different choice of initial condi- 
tions. They are even different for uncorrelated and weakly correlated initial conditions with 
the same initial concentration. The corrections are smaller for the interpolating sequence of 
particles and holes than for the random distribution in the case of open boundary conditions. 

4 Conclusions 

In this paper, we presented the numerical examination of the finite-size scaling behaviour 
of the coagulation model. 

We have shown on the one hand the universality of the scaling function with the help of 
Monte Carlo simulations. On the other hand, we were able to determine critical exponents 
already from lattices of 10 sites using finite-lattice extrapolations. The test of these two 
numerical methods (Monte-Carlo simulations and a finite-lattice extrapolations) revealed that 
they can be successfully applied to characterize the finite-size scaling behaviour of reaction- 
diffusion processes. 

Let us discuss our results in detail. 

• Universality of the scaling function and the corrections 
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As a first application of the Monte Carlo simulations, we examined the scaling function 
for periodic and open boundary conditions for different tuning of the rates in order to 
answer the question whether the scaling function is universal (i.e. independent of the 
details of the model): we varied the parameter A' that reflects the difference between 
the diffusion rate (D) and the coagulation rate (c). The scaling function was found to 
be universal. 

As far as the corrections to the scaling function are concerned, we obtained a leading 
correction term of the order in the case of A' 7^ 0. The corresponding correction 
function (the coefficient of is different for different values of A'. However, it has 
always the same sign as A' (cf. Eq. ( [2.1QQ ). Especially, for A' = it is zero as well. 
This explains the fact that we obtained leading corrections of the order -p- in the case 
A' = which was treated analytically in paper I and II. 

• Influence of the initial conditions 

We examined the influence of weakly correlated initial states both for open and for 
periodic boundary conditions. A configuration with small clusters (AA$$AA . . .) was 
used as initial state. We have shown that the scaling function is not influenced by the 
initial conditions. So we observe again the phenomenon of self-organisation. The system 
develops according to its intrinsic dynamics and is independent of the initial conditions 
for large times and large lattices. 

• Test of numerical methods 

In order to test the accuracy of Monte-Carlo simulations, we compared the results ob- 
tained by simulation to the data that are analytically known. This comparison revealed 
that the Monte-Carlo method provides reliable results that approximate very well the 
exact expressions. So this method can successfully be used to investigate the finite-size 
scaling behaviour of chemical models. 

As far as finite lattice extrapolations are concerned, we have demonstrated that the 
algorithm of Bulirsch and Stoer || is a useful method to determine the scaling and the 
correction exponents and -functions already from very small lattices. The only condition 
one has to impose is that the input values (for small lattices) have to be known with 
a sufficiently high precision. The accuracy of the extrapolated values depends on the 
number of input values. However, in general it is sufficient to take the values for 8—10 
sites as input values to determine the exponents with a precision of at least 10~ 3 . So this 
method is very suitable for the analysis of the finite-size scaling regime for models that 
cannot be treated analytically. In this case a numerical diagonalisation is still possible 
for small lattices and can be used to generate the input values for the BST-extrapolation. 

• Best choice of the parameters for the numerical examination 

The analysis of the coagulation model allows to decide which is the best choice of the 
parameters (tuning of the rates, initial conditions, boundary conditions and value of the 
scaling variable z) for the numerical determination of the finite-size scaling data. 

1. Tuning of the rates 
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As far as this model is concerned, the scaling function can best be determined if 
the rates are tuned in a way that A' = 0. Then the leading correction term is 
of the order L~ 2 , i.e. small for large L; furthermore the correction function (the 
coefficient of L~ 2 ) is smallest for A' = 0. 

2. Choice of initial conditions 

As far as the choice of the initial conditions is concerned, the corrections become 
smallest for large values of the initial occupation probability p as can be seen 
in Fig. 5. The effect of correlated initial conditions (an interpolating sequence 
of particles and holes) cannot be easily understood and requires a more detailed 
analysis: for periodic boundary conditions, the corrections become smaller for 
uncorrelated initial conditions (when the initial concentration is kept the same) 
while for open boundary conditions, the cluster configuration gives rise to smaller 
correction terms (cf. Tabs. [X] and [XII] ). So in general the most promising choice 



of initial conditions for the determination of the scaling function and the scaling 
exponent is the full lattice. 

3. Choice of boundary conditions 

The discussion which boundary conditions are more convenient for a numerical 
examination of the scaling and the correction exponent (the scaling functions are 
different for different boundary conditions) is very difficult and depends on the 
other parameters as well. If the initial occupation probability p is smaller than 
1 and A' = the corrections are smaller for periodic boundaries. For A' ^ 
and p = 1 however, the corrections are smaller for open boundary conditions 
(cf. Tabs. |V| and |V11J| ). These cases are interesting for the determination of the 



correction exponent. Taking p = 1 and A' = (the best choice of A' and p for 
the determination of the scaling exponent) the corrections are almost identical for 
different boundary conditions. In Tabs. [V| and [VIII| values for z between 0.01 and 



0.1 are given where the corrections are exactly the same. For larger values of z, 
however, they become slightly different (cf. paper II). Taking into account that 
periodic boundaries usually allow a Fourier transformation of the equations and 
therefore a reduction of the number of degrees of freedom of the system, it is more 
convenient to work with periodic boundary conditions here. 

4. Value of the scaling variable z 

For the determination of the scaling exponent, one has to choose a value of z where 
the corrections are small, i.e. a value around z = 0.6 should be convenient. Since 
the zero of the correction function occurs at a smaller value of z for periodic bound- 
ary conditions, these boundary conditions better allow a Monte Carlo simulation: 
the smaller the value of z is, the smaller the simulated time can be made and 
the more precise the results of the simulation are. For the determination of the 
correction exponent the situation is quite different. Here the correction function 
should be large in order to allow a numerical fit of the simulated data. So one 
has to choose small values of z for a Monte Carlo simulation (z < 0.02). The 
BST-algorithm, however, converges better for large values of z (z > 0.2). 

Whether these observations hold for more complicated chemical systems still has to be 
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investigated. 

Summing up, we can say that the numerical treatment presented here completes the 
picture of the finite-size scaling behaviour of the coagulation model. Both, analytical and 
numerical methods can be succesfully applied to chemical models. The results obtained in 
our work now open the way towards the investigation of more complex chemical systems. 
The methods tested here provide a powerful tool to gain deeper insight into non-equilibrium 
physics. 
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At 


L 


tmax 


Maximum error in %. 


5- 10- 3 


6 


0.9 


0.36 


5-10- 2 


200 


10 3 


0.42 


5- 10- 1 


1000 


25 • 10 3 


0.44 



Table I: Maximum error of the simulated concentration for periodic boundaries and A' = 





Scalin; 


5 function F (z) 


Corrections Fi(z, 0) 


z 


anal. 


M.C. 


anal. 


M.C. 


0.01 


7.979 


7.839 ± 0.005 


-99.8 


-103.2 ±0.3 


0.02 


5.642 


5.680 ±0.004 


-35.3 


-45.0 ±0.3 


0.03 


4.607 


4.635 ± 0.004 


-19.2 


-24.3 ±0.3 


0.04 


3.989 


4.006 ±0.003 


-12.5 


-15.0 ±0.3 


0.05 


3.568 


3.577 ±0.003 


-8.9 


-10.1 ±0.3 


0.06 


3.257 


3.264 ±0.003 


-6.8 


-7.4 ±0.2 


0.07 


3.016 


3.022 ±0.003 


-5.4 


-5.8 ±0.2 


0.08 


2.821 


2.825 ±0.003 


-4.4 


-4.7 ±0.2 


0.09 


2.660 


2.661 ±0.003 


-3.7 


-3.7 ±0.2 


0.10 


2.523 


2.525 ±0.003 


-3.2 


-3.2 ±0.2 



Table II: Monte Carlo data and exact results for periodic boundaries and A' = 



z 


A' = 0.75 


A' = 0.5 


A' = 


A' = -1 


A' = -2 


0.01 


7.95 ±0.06 


7.87 ±0.03 


7.98 


7.97 ±0.02 


7.88 ±0.03 


0.02 


5.43 ± 0.04 


5.58 ±0.02 


5.64 


5.62 ±0.02 


5.57 ±0.02 


0.03 


4.50 ±0.04 


4.55 ±0.02 


4.61 


4.58 ±0.02 


4.60 ±0.02 


0.04 


3.91 ±0.04 


3.94 ±0.02 


3.99 


3.95 ±0.02 


3.96 ±0.02 


0.05 


3.50 ±0.03 


3.55 ±0.02 


3.57 


3.54 ±0.02 


3.56 ±0.02 


0.06 


3.20 ±0.03 


3.25 ±0.02 


3.26 


3.23 ±0.02 


3.24 ±0.02 


0.07 


2.97 ±0.03 


3.00 ±0.02 


3.02 


2.98 ±0.02 


3.00 ±0.02 


0.08 


2.78 ±0.03 


2.80 ±0.02 


2.82 


2.80 ±0.02 


2.79 ±0.02 


0.09 


2.64 ±0.03 


2.63 ±0.02 


2.66 


2.63 ±0.01 


2.63 ±0.02 


0.10 


2.49 ±0.03 


2.51 ±0.02 


2.52 


2.51 ±0.01 


2.50 ±0.02 



Table III: Scaling function F (z) for periodic boundaries 
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z 


A' — n 7^ 

L\ — U. ( 


A' n x. 


A' — n 

L\ — U 


A' — 1 


A' — 9 


n m 
U.Ul 


4 ( 1 ± 




n 
U 


KCi 7-1-0/1 

— oy. 1 ± z.4 


— ( U ± Z.U 




251 ± 5 


75 6 ± 2 6 

( '.J . \f _i_ _ . \j 


n 


—28 8 ± 2 


—35 ± 1 8 


0.03 


160 ±4 


51.2 ±2.3 





-18.3 ± 1.9 


-26 ± 1.6 


0.04 


118 ±4 


38.2 ±2.1 





-11.7 ± 1.8 


-18 ± 1.5 


0.05 


93 ±4 


28.8 ±2.0 





-9.9 ± 1.7 


-16 ± 1.5 


0.06 


76 ±3 


22.6 ± 1.8 





-7.9 ± 1.6 


-13 ± 1.4 


0.07 


65 ±3 


20.0 ± 1.7 





-5.7 ± 1.5 


-11 ± 1.4 


0.08 


57 ±3 


19.0 ± 1.7 





-5.4 ± 1.5 


-8 ±1.3 


0.09 


49 ±3 


18.0 ± 1.6 





-4.2 ± 1.5 


-7± 1.3 


0.10 


45 ±3 


14.3 ± 1.6 





-4.8 ± 1.4 


-6 ±1.3 



Table IV: First order correction function F ± (z, A') for periodic boundaries 



z 


A' = 0.75 


A' = 0.5 


A' = 


A' = -1 


A' = -2 


0.01 


-8212 ± 156 


-2137 ± 74 


-99.8 


448 ± 53 


418 ± 38 


0.02 


-3197 ± 124 


-766 ± 60 


-35.3 


147 ± 46 


148 ± 33 


0.03 


-1516 ± 108 


-438 ± 52 


-19.2 


56 ±42 


124 ± 30 


0.04 


-922 ± 97 


-289 ± 48 


-12.5 


-16 ±39 


61 ±28 


0.05 


-634 ± 90 


-168 ±45 


-8.9 


-9 ±37 


72 ±27 


0.06 


-434 ± 84 


-97 ±42 


-6.8 


-20 ± 36 


50 ±26 


0.07 


-332 ± 80 


-84 ± 40 


-5.4 


-43 ± 34 


29 ±25 


0.08 


-273 ± 76 


-97 ±39 


-4.4 


-35 ± 33 


4 ±24 


0.09 


-177 ±73 


-121 ±38 


-3.7 


-46 ± 33 


-4 ±23 


0.10 


-183 ±71 


-68 ± 36 


-3.2 


-15 ±32 


-7 ±23 



Table V: Second order correction function F 2 (z,A') for periodic boundaries 



z 


A' = 0.5 


A' = 


A' = -1 


0.01 


8.33 ±0.04 


8.34 


8.27 ±0.05 


0.02 


5.98 ±0.03 


6.01 


6.00 ±0.04 


0.03 


4.92 ±0.03 


4.97 


4.99 ±0.04 


0.04 


4.28 ±0.02 


4.35 


4.41 ±0.04 


0.05 


3.89 ±0.02 


3.93 


3.96 ±0.03 


0.06 


3.60 ±0.02 


3.62 


3.65 ±0.03 


0.07 


3.36 ±0.02 


3.38 


3.40 ±0.03 


0.08 


3.15 ±0.02 


3.18 


3.18 ±0.03 


0.09 


3.01 ±0.02 


3.02 


3.01 ±0.03 


0.10 


2.88 ±0.02 


2.89 


2.85 ±0.03 



Table VI: Scaling function F Q (z) for open boundaries 
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z 


A' — n ^ 

Z-i — u.o 


A' — n 


A' — 1 


n m 


132 ± 3 


n 


—58 + 4 

KJ _L_ t: 


0.02 


65 ±2 





-33 ±4 


0.03 


44 ±2 





-25 ±3 


0.04 


36 ±2 





-24 ±3 


0.05 


27 ±2 





-17 ±3 


0.06 


21 ±2 





-15 ±3 


0.07 


18 ±2 





-13 ±3 


0.08 


16 ±2 





-9 ±3 


0.09 


13 ±2 





-7 ±3 


0.10 


11 ±2 





-4 ±3 



Table VII: First order correction function Fi(z,A') for open boundaries 



z 


A' = 0.5 


A' = 


A' = -1 


0.01 


-1845 ±56 


-99.7 


355 ± 97 


0.02 


-648 ± 45 


-35.3 


180 ± 83 


0.03 


-386 ± 40 


-19.2 


173 ± 76 


0.04 


-308 ± 36 


-12.5 


215 ±71 


0.05 


-207 ±34 


-8.9 


123 ± 68 


0.06 


-144 ±32 


-6.8 


125 ± 66 


0.07 


-112 ±31 


-5.4 


97 ±63 


0.08 


-113 ±30 


-4.4 


29 ±62 


0.09 


-71 ±29 


-3.7 


7 ±60 


0.10 


-55 ± 28 


-3.2 


-39 ± 59 



Table VIII: Second order correction function F 2 (z,A') for open boundaries 



A' 


a(A', 1) 


-2.00 


-0.512 ±0.060 


-1.00 


-0.414 ±0.083 


0.50 


0.894 ±0.070 


0.75 


2.831 ±0.130 



Table IX: First correction term a(A', 1) of the large-time expansion for periodic boundaries 
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Scaling 


function Fq(z) 


Corrections F 2 (z,0) 




p = 0.5 


AAM . . . 


p = 0.5 


AAM . . . 


z 


anal. 


M.C. 


anal. 


M.C. 


0.01 


7.979 


7.967 ±0.008 


-212.5 


-836.7 ± 12.6 


0.02 


5.642 


5.623 ±0.007 


-106.3 


-304.8 ± 10.7 


0.03 


4.607 


4.607 ±0.006 


-70.8 


-164.4 ±9.7 


0.04 


3.989 


3.988 ±0.006 


-53.1 


-109.8 ±9.0 


0.05 


3.570 


3.556 ±0.006 


-42.5 


-77.2 ±8.6 


0.06 


3.257 


3.245 ± 0.005 


-35.4 


-61.7 ±8.2 


0.07 


3.016 


3.006 ±0.005 


-30.4 


-48.8 ± 7.9 


0.08 


2.821 


2.808 ±0.005 


-26.6 


-40.6 ± 7.6 


0.09 


2.660 


2.647 ±0.005 


-23.6 


-27.7 ± 7.4 


0.10 


2.523 


2.517 ±0.005 


-21.3 


-17.6 ±7.2 



Table X: Influence of the initial conditions on the scaling properties for periodic boundaries 
and A' = 





l + x 


Scaling function F (z) 


-2-y 


Corrections F 2 (z) 


z 


extr. 


extr. 


analytical 


extr. 


extr. 


analytical 


0.1 


0.004547 


2.5210 


2.5231 


-0.134 


-2.776 


-3.154 


0.2 


0.002986 


1.784283 


1.784286 


-0.011 


-1.1011 


-1.1017 


0.3 


0.0000353 


1.46043902 


1.46043901 


-0.005 


-0.51697 


-0.51697 


0.4 


0.0000003 


1.27857491 


1.27856700 


0.005 


-0.22442 


-0.22246 


0.5 


0.0000477 


1.169713390 


1.169713392 


0.026 


-0.07196 


-0.07196 


0.6 


-0.000002 


1.103560880 


1.103560906 


0.145 


-0.00291 


-0.00291 


0.7 


0.0000304 


1.0632169524 


1.0632169522 


-0.003 


0.023694 


0.023694 


0.8 


0.0000726 


1.0385928840 


1.0385928831 


0.005 


0.030097 


0.030097 


0.9 


0.0000522 


1.02356074810 


1.02356074816 


-0.002 


0.0279317 


0.0279317 


1.0 


0.0000988 


0.01438377205 


0.01438377206 


0.006 


0.0228893 


0.0228893 



Table XI: Extrapolants from finite lattices (L — 2, ... ,9) for periodic boundary conditions 
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Scaling function F (z) 


Corrections F 2 (z) 


z 


p = 0.5 


AAM . . . 


p = 0.5 


AAM . . . 




anal. 


extrapol. 


anal. 


extrapol. 


0.01 


8.342226 


8.344116 


-1950.15 


-1050.17 


0.02 


6.005276 


6.005276 


-726.78 


-381.02 


0.04 


4.352803 


4.352803 


-275.60 


-144.03 


0.07 


3.379100 


3.379100 


-127.93 


-66.66 


0.10 


2.886513 


2.886513 


-79.08 


-41.12 


0.20 


2.147403 


2.147403 


-31.71 


-16.41 


0.30 


1.818125 


1.818125 


-18.94 


-9.79 


0.40 


1.617195 


1.617195 


-13.21 


-6.81 


0.50 


1.475872 


1.475872 


-9.83 


-5.04 


0.60 


1.369969 


1.369969 


-7.50 


-3.82 


0.70 


1.288534 


1.288534 


-5.77 


-2.91 


0.80 


1.225287 


1.225287 


-4.45 


-2.23 


0.90 


1.175981 


1.175981 


-3.44 


-1.70 


1.00 


1.137489 


1.137489 


-2.66 


-1.30 



Table XII: Influence of the initial conditions on the scaling properties for open boundaries 
and A' = 
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